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We explore in detail the possibility of generating a pair-coherent state in the non-degenerate 
parametric oscillator when decoherence is included. Such states are predicted in the transient regime 
in parametric oscillation where the pump mode is adiabatically eliminated. Two specific signatures 
are examined to indicate whether the state of interest has been generated, the Schrodinger cat 
state — like signatures, and the fidelity. Solutions in a transient regime reveal interference fringes 
which are indicative of the formation of a Schrodinger cat state. The fidelity indicates the purity of 
our prepared state compared to the ideal pair-coherent state. 



I. INTRODUCTION 

Quantum mechanics is well known for two fundamen- 
tal quantities, the superposition of states and entangle- 
ment. Both have been studied extensively and much rich 
physics has been observed. Generally however systems 
exhibit one but not generally both of these phenomena. 
Recently however the idealised state produced in non- 
degenerate parametric oscillation when the pump mode 
is adiabatically eliminated has been shown to give the 
possibility of observing both Schrodinger cat states |^] 
and producing correlated photon pairs which could be 
used to perform tests of quantum mechanics versus lo- 
cal realism such as the EPR paradox ||] and the Bell 
inequality |,|. 

There is still much interest in the possibility of generat- 
ing experimentally a Schrodinger cat state. Such a state 
is defined as a quantum superposition of two macroscopi- 
cally distinct states It was Schrodinger's concern that 
quantum mechanics does not fundamentally prohibit the 
existence of such states, which seemingly defy physical 
reality. Work by Krippner and Reid predicts that a 
Schrodinger cat state may be produced in the signal field 
of non-degenerate parametric oscillation, in a transient 
regime. The theoretical model presented here incorpo- 
rates the effect of linear signal losses, which tend to op- 
pose the formation of the superposition state. 

Recent work by Gilchrist et. al. and Munro ||] 
showed how in a certain narrow regime, the pair-coherent 
state gives quantum mechanical predictions that are in 
disagreement with those of local hidden variable theo- 
ries for a situation involving continuous quadrature phase 
amplitude measurements. This test could be achieved by 
binning the continuous position and momentum informa- 
tion into two categories and using the binary results in 
the strong Clauser Horne Bell inequality test |Q. The 
predicted violation was small (less than 2%). Such re- 
sults were highly idealised, and assumed the preparation 
of a pair-coherent state to begin with. 

In this paper we will investigate the possibility of gen- 



erating a pair-coherent state in the non-degenerate para- 
metric oscillator with an adiabatically eliminated pump 
mode when decoherence is included. Two specific signa- 
tures of the state will be considered, interference fringes 
which are indicative of the formation of a Schrodinger 
cat state and the fidelity. In particular the fidelity will 
provide an indication of the parameter regime required 
to perform a loophole free test of quantum mechanics. 

It should be noted that the fidelity is a rather ab- 
stract quantity and in practice would require the en- 
tire state to be reconstructed in order to be measured. 
The Schrodinger cat state signatures however are more 
operationally accessible though given the sensitivity of 
Schrodinger cat state itself to decoherence Q we would 
expect similar sensitivity for the pair-coherent state and 
this is indeed the case. 



II. THE PAIR-COHERENT STATE 

The following two-mode entangled quantum superpo- 
sition state, 

|circle),„ = N r e-''™^>oe'^")a|roe-''^)fc* (1) 
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is known as the was pair-coherent (or "circle" ) state and 
was originally discussed by Agarwal [p|jlo[| and Reid and 
Krippner In equation (|l|) |. . .)a and |. . .);, represent 
coherent states in the modes a and h, where the a and 
b are the usual boson operators. A/" is a normalisation 
coefficient and the amplitude of the coherent state, m 
is the photon number difference between the signal and 
idler modes. This state is actually a continuous super- 
position of coherent states in a circle (hence the notation 
I circle)). For our purposes in this paper we shall concen- 
trate on the m = (equal photon number in each mode) 
situation in equation (|l]) which when normalised is given 
by 
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where 
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Here Jq is a zeroth order modified Bessel function. Such 
a state can also be written in terms of correlated photon 
number pairs of the form 



|vl/) = ^c„|n)|n) 
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where 



J2n 
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Such a state should not be confused with the state 
produced by the non-degenerate parametric amplifier 
(NDPA) which can also be written in the form of (Q) 
but now with the c„ coefficients given by 



tanh" [xer] 
cosh [x^t] 



(6) 



Here e represents the field amplitude of a non-depleting 
classical pump, x is proportional to the susceptibility of 
the medium and r is the time that the modes spend in 
the crystal. 

In this paper we are interested in exploring the gener- 
ation of the pair-coherent state. 



III. THE NON-DEGENERATE PARAMETRIC 
OSCILLATOR 

It has been suggested by Reid and Krippner that the 
NDPO transiently generates a state of the above form, in 
the limit of very large parametric nonlinearity and high- 
er cavities [Q. The non-degenerate parametric oscillator 
with linear damping in the signal and idler modes and 
the pump mode pumped by a classical field may be rep- 
resented by the Hamiltonian 



H — Hi + Hp + -ffincv 
Hi = ihK{a3a2a\ — a^^a2ai) 
H„ 



(7) 
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Here are boson operators for the cavity modes at fre- 
quencies oji, where ujj, ^ uji + uj2- The mode 03 is driven 
by a resonant external driving field with amplitude pro- 
portional to e and is known as the pump-mode. Modes ai 
and a2 are the signal- and idler-modes. The loss of pho- 
tons through the cavity mirrors is modeled by the Hamil- 
tonian term TJjrrcv, which denotes a coupling of the cavity 



modes to the zero temperature reservoir modes (symbol- 
ised by Ti) external to the cavity. We will denote the 
cavity decay rates for the modes ai by 7^. 

Following standard techniques it is easy to derive a 
master equation of the form 



p^^[Hi + Hp,p\ 
in 
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In the limit where the pump mode is heavily damped 
compared to the other modes (73 ^ 72,71) the pump 
variables can be eliminated. This is equivalent to study- 
ing the following model Hamiltonian: 

6Ai /ft 
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The presence of the two-photon damping term is the fun- 
damental difference between the NDPO in this limit and 
the NDPA. 

To aid our discussion of parameters below, we will 
briefly examine realistic parameter values for the non- 
degenerate parametric oscillators containing the com- 
monly used crystals, silver gallium selinide (AgGaSe2) 
and potassium titanyl phosphate (KTP). In Table (|) 
are shown some typical values. We easily observe that 
the nonlinear coupling constant is much weaker than 
the damping constant and hence our scaled parameter 

= K^/773 that we will introduce shortly will be very 
small. 



CRYSTAL 




7(s~ 




k/7 


AgGaSe2 


4.4 X 10* 


7.5 X 




5.9 X 10"^ 


KTP 


7.6 X 10^ 


7.5 X 


10** 


1 X 10"^ 



TABLE I. Table of realistic values for the nonlinear cou- 
pling constant k and the damping constant 7 for two types of 
parametric crystal AgGaSe2 and KTP. 
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IV. THE ADIABATICALLY ELIMINATED 
MASTER EQUATION 



V. ENTANGLEMENT 



The Hamiltonian above corresponds to the follow- 
ing master equation (equation with the pump mode 
adiabatically eliminated) 



CIT 



—g'^(2a\a\paia2 — a\a\aia2P — pa| 020102) 
2 

+ ^(2a]pOj - o]ajp - pa^^aj) (10) 



where we have introduced the following scaled variables 

A = and q = The time has been scaled such 

737 ^ 737 

that T — and we have assumed that the signal and 
idler decay constant 71 and 72 are in fact equal to 7. 

The master equation can be solved numerically by pro- 
jectin g th e master equation onto an infinite number state 
basis 11 1 . Expanding the density matrix in the number 
state basis as follows 

{ni\{n2\p\mi)\m2) , (11) 
we may express the time evolution of the system as 
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where this super matrix C, 



nin2 ;mim2 



is given by 



nin2;mim2 



rni + l,mi rn2 + l,m2 
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As we mentioned in the introduction, this pair- 
coherent state has the property that it contains sufficient 
entanglement to violate a Bell inequality. Gilchrist et. al. 
showed that the pair-coherent state specified by (P the- 
oretically violated a Bell inequality. To be more explicit, 
they showed how using highly efficient quadrature phase 
homodyne measurements, the Clauser Horne strong Bell 
inequality could be tested in an all optical regime. While 
the violation may be small the highly efficient detection 
means that provided the extremely ideal state could be 
generated, a significant test could be done. 

There are a number of measures to determine the pu- 
rity of the produced state. The measure we will use here 
is the fidelity. The fidelity may be defined as 



F = |(circle|output)|^ 



(15) 



in terms of pure states. In terms of the density operator 
p of the output state, we represent the fidelity as 



F = Tr 



1/2 
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Figure |i| shows the result of calculating the fidelity 
against an ideal pair-coherent state for X/g^ = 1.12 and 
two different values of It shows a maximum fidelity 
of around 80% for = 300. For larger g we can get a 
larger fidelity but the transient period over which this is 
available is significantly shorter. Given the small size of 
the Bell inequality violation, and the narrow parameter 
regime over which it occurs in reference it is likely that 
higher nonlinearities would be required for the NDPO to 
produce the state sufficient to violate the Bell inequality 
in that scheme. 

Gilchrist: Figure 1 
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Here 



1 if i = n and j = m 
otherwise 



(14) 



To allow for numerical calculations, one must put a fi- 
nite limit on the number of Fock states used in the basis 
in ( p^ . Care must be taken to ensure that the truncation 
of the number state basis is done correctly so the popu- 
lation of the higher order states is small. In practice we 
found nmax = 20 sufficient for most of the calculations. 




0.10 



FIG. 1. Plot of the fidelity of the generated state compared 



with the pair-coherent state for \/g^ 
giving a maximum of F ~ 0.78, (ii) 
mum of F ~ 0.76. 



.12 and (i) = 300 
= 10 giving a maxi- 
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VI. SCHRODINGER CAT STATES 

Given a pair-coherent state, previous work by Krippner 
and Reid 0| has predicted that in the hmit of very large 
ro and with a conditional measurement of one mode, 
a Schrodinger cat state may be produced in the other 
mode. Observing these states then constitutes an indi- 
rect signature of the presence of a pair-coherent state as 
well as being of interest in its own right. In this sec- 
tion we examine the formation of these states when the 
parameters are not so extreme and with the presence of 
damping. 

Let us now suppose that one measures the quadrature 
phase amplitude defined by 



IV2 



where 9i is the phase of the local oscillator for the zth 
mode. The two measured quadratures Xq and are 
non commuting observables. We assume that the mea- 
surement of the quadrature Xg^ gives the result xg^ . 

It is then possible to construct the joint probability 
distribution of obtaining a result xg-^ for the first mode 
and xg^ for the second. This probability is expressed as 



Pet,e2 {.xi,X2) = {xi\{x2\p\x2)\xi 



(18) 



In figure g we plot Pq.o (a^i = z,X2 = 0) versus z and 
[xi = p,X2 = 0) versus p. We clearly observe the 



P. 



7r/2,0 

interference fringes and twin peaks that characterise the 
Schrodinger cat state \iX/ g^) + | — iX/g^). 




FIG. 2. Plot of position (Xq) (a) and momentum 
(b) probability distribution for the signal mode conditioned on 
a quadrature position measurement on the idler mode record- 



(17) ing 0. Here we have \/g^ — 1.5. 



For comparison, the Schrodinger cat state signatures 
are also shown for the pair-coherent state required by 
Gilchrist et. al. (A/.g^ — 1.12). As can be expected 
we do not get very distinct interference fringes and peak- 
resolution. Interestingly, clear cat signatures require only 
slightly larger values of \/g^. 



Here we have abbreviated xg-^ by xi and xg^ by X2- 
Our state of interest can then be detected by observ- 
ing interference fringes in the probability distribution 
-Pir/2,0 (a^i = z,X2 = 0) for the quadrature phase ampli- 
tude measurement performed on the signal mode, condi- 
tioned on the idler mode result X2 = for O2 = 0. The 
observation of interference fringes present in one of the 
quadrature measurements (in conjunction with the obser- 
vation of twin isolated peaks in the conjugate quadrature 
phase amplitude) are indicative of Schrodinger cat states, 
where Xq is analogous to the role of position and 
of momentum. 

In terms of our number state basis expansion for the 
density matrix, this joint probability distribution can be 
written as 

00 2 
Pg^fi^{xi,X2) = ^ ]^(xeJn)(TO|xeJ (19) 





FIG. 3. Plot of position (Xo) (a) and momentum (Jf^/2) 
(b) probability distribution for the signal mode conditioned on 
a quadrature position measurement on the idler mode record- 
ing 0. Here we have \/g^ — 1.12. 



B. Numerical Simulations 



The probability distribution is then written as 



where {xg^\n) is given by 
{xg\n) = 



Pg 



{X1.,X2) 



E 
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im2 Yli^Bi \n) {in\xg^) (21) 
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Hn (xg) 
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In our calculations the effect on increasing the number of 
basis states by one produced an error of less than 0.001 
percent. 

The results of our calculations are shown in the Fig- 
ures Figure ^ plots the position (Xq) and momentum 
(^ir/2) probability distributions versus time. The inter- 
ference fringes in the momentum probability distribution 
In the absence of damping the ideal pair-coherent state combined with the twin peaks in the position distribution 
is given by (||) with the c„ coefficients specified by (||). reveal the Schrodinger cat state-like nature. 



where the units have been chosen such that h = lo = c 
1 and Hn (xg) is the Hermite polynomial. 

A. The ideal situation 
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Gilchrist: Figure 5a 
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FIG. 4. Plot of (a) the evolution of the position probability 
distribution P{z), and (b) the momentum probability distri- 
bution P(p). Here = 300 with X/g'^ = 1.5. 
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Gilchrist: Figure 5b 
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FIG. 5. Plot of (a) the position probability distribution 
P{z) and (b) the momentum probability distribution P(p); 
for (i) g'^ = 300 with t = 0.014r, (ii) g"^ = 10 with t = 0.059r, 
(iii) g^ — 3 with t — 0.19r. Here X/g^ = 1.5, and the time 
chosen was informally optimised. 



The formation of fringes with the evolution of the sig- 
nal field from the vacuum state is clearly evident in Fig- 
ure ^. As the oscillator evolves further the fringes are 
washed out. The \iX/g'^) — \ — iX/g'^) state, which is gen- 
erated from |«A/(jf^) -I- 1 — iX/g"^) with the loss of a cavity 
photon, contributes more significantly as time increases, 
and the fringes are lost in this case after only O.lr. 



In order to establish the orders of g required to obtain 
a clear fringe pattern, the P{z) and P{p) distributions 
are shown in figure |^ for a range of g with X/ g^ — 1.5. 
For g^ greater than or of the order of 10, interference 
fringes become apparent in the transient evolution of the 
oscillator. The fringes (for fixed X/g'^) become more pro- 
nounced as g increases. This is consistent with the earlier 
analytical conclusions, which were based on calculations 
performed in the large g limit where the strength of the 
two-photon nonlinearity is much greater than the single- 
photon cavity loss rate. 



As a final appraisal, we can test the fidelity of the cat 
state against an ideal Schrodinger cat state of the form 
\iX/ g"^) ~\- 1 — iX/g^) and this is shown in figure ^ for var- 
ious times. Note that it is perfectly possible to generate 
a Schrodinger cat state not of this form but which still 
constitutes a superposition of two macroscopically dis- 
tinct states, hence the fidelity in this case is an indication 
of the purity of the underlying pair-coherent state rather 
than an indication of a good Schrodinger cat state. 
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Gilchrist: Figure 6 




0.0 L , , \ , , , \ , , , I , , , \ , , ^ 

0.00 0.02 0.04 0.06 0.08 0.10 

T 

FIG. 6. Plot of the fidelity of the generated cat state for 
(i) = 300 giving a maximum of F ~ 0.94 and (ii) — 10 
giving a maximum of ~ 0.78. Here X/g^ = 1.5. The result 
is conditioned on a quadrature position measurement on the 
idler mode recording 

In this section we have established how the non degen- 
erate parametric oscillator operating in an adiabatically 
eliminated pump mode regime will produce Schrodinger 
cat states in the transient evolution. 



VII. CONCLUSIONS 

In this paper we have examined two signatures of the 
formation of pair coherent states in the NDPO in the 
presence of linear damping. The most direct indication 
is to look at the fidelity and even with large nonlincari- 
ties we found the fidelity to be poor although higher non- 
linearity improves the fidelity. Certainly the acquired fi- 
delity would indicate that the generated state would be a 
poor candidate for the Bell inequality test of Gilchrist et. 
al. even for nonlinearities as high as given by ~ 300. 



A more indirect measurement of the purity of the state 
is to look for the formation of Schodinger cat states which 
are predicted upon a conditioned measurement on one 
mode given a large value for \/g^ . The formation of these 
states not only gives an indication of the presence of a 
pair-coherent state but are of interest in their own right. 
Here we predict that for nonlinearities characterised by 

~ 300 formation of clear Schodinger cat states is pos- 
sible for only X/ g^ ^ 1.5. Though, again, this parameter 
regime is difficult to produce experimentally. 
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